A note on digitized angles 
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Abstract. We study the configurations of pixels that occur when two digitized 
straight lines meet each other. 

About ten years ago I was supervising the Ph.D. thesis of Chris Van Wyk [4], which 
introduced the IDEAL language for describing pictures [5] . Two of his example illustrations 
showed arrows constructed from straight lines something like this: 



When I looked at them, I was sure that there must be a bug either in IDEAL or in the 
TROFF processor that typeset the IDEAL output, because the long shafts of the arrows did 
not properly bisect the angle made by the two short lines of the arrowheads. The shafts 
seemed to be drawn one pixel too high or too low. Chris spent many hours together with 
Brian Kernighan trying to find out what was wrong, but no errors could be pinned down. 
Eventually his thesis was printed on a high-resolution phototypesetter, and the problem 
became much less noticeable than it had been on the laser-printed proofs. There still was 
a glitch, but I decided not to hold up Chris's graduation for the sake of a misplaced pixel. 

I remembered this incident at the end of 1983, when I was getting ready to write a 
new version of the METRFONT system for digital art [3]. I didn't want my system to have 
such a flaw. But to my surprise, I learned that the problem is actually unavoidable in 
raster output: It is almost impossible to bisect a digitized angle exactly, except in very 
special circumstances. The two "halves" of the angle will necessarily appear somewhat 
different from each other, unless the resolution is quite high. Therefore Van Wyk (and 
Kernighan) were vindicated. Similar problems are bound to occur in MacDraw and in any 
other drawing package. 

For example, one of the things I noticed was the following curious fact. Consider the 
45° angle that is made when a straight line segment of slope 2 comes up to a point (xq, yo) 
and then goes down along another line of slope —3: 




( Xo + 2t,y -t)/ \ {xQ + u ^ yo _ 3u) 



If we digitize this angular path, the result will take one of five different shapes, depending on 
the value of the intersection point (xq, yo), whose coordinates are not necessarily integers. 
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The possibilities are 




Pi 






Pa 




If this 45° angle provides the left half of an arrowhead, the right half of the arrowhead will 
be a 45° angle made by a line of slope —3 meeting a line of slope —1/2, 

Od»yi) 

\ ^^(a;i + 2it, yi - u) 
(xi +t,y 1 - 3t)\ 

For this angle there are, similarly, five possibilities after digitization, namely 




Qi 






Qa 




To complete the arrowhead, we should match the left angle Pi with an appropriate Qj. 
But none of the Q's has the same shape as any of the P's. And this is the point: Human 
eyes tend to judge the magnitude of an angle by its appearance at the tip. By this 
criterion, some of these angles appear to be quite a bit larger than others (except at high 
resolutions). Hence it is not surprising that a correctly drawn angle of type P would appear 
to be unequal to a correctly drawn angle of type Q, even though both angles would really 
be 45° when drawn with infinite resolution. (The patterns of white pixels, not black pixels, 
are the source of the inconsistency.) Here, for example, are four quite properly digitized 
arrows with shafts of increasing thickness: 



xH* list; 



xxxxxx 



xxx 
XX 



xxx 
xxx 
xxx 

XXX 

xxx 
xxx 
xxx 



Hi!** 
* *!!! 

xxx* 

XX 



xxxxxx 
*xxxx 

XX 



We might also want to know the probability that the digitized shape will be of a 
particular type when the corner point (xo,yo) is chosen at random in the plane. Is 
one of the patterns more likely to occur than the others? The answer is no, when we use 
the most natural method of digitization; each Pj. will be obtained with probability 1/5. 
Similarly, each of the five shapes Qk turns out to be equally likely, as (xi, y±) varies. 

The main purpose of this note is to prove that the facts just stated are special cases 
of a general phenomenon: 

Theorem. When a line of slope a/b meets a line of slope c/d at a point (xo,yo), the 
number of different digital shapes it can produce as (xq, yo) varies is \ad — bc\. Moreover, 
each of these shapes is equally likely to occur, if (xq, yo) is chosen uniformly in the plane. 

We assume that a/b and c/d are rational numbers in lowest terms. Two digitized shapes 
are considered to be equal if they are identical after translation; rotation and reflection are 
not allowed. 
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Before we can prove the theorem, we need to define exactly what it means to digitize 
a curve. For this, we follow the general idea explained, for example, in [3, Chapter 24]. 
We consider the plane to be tiled with pixels, which are the unit squares whose corners 
have integer coordinates. Our goal is to modify a given curve so that it travels entirely on 
the boundaries between pixels. If the curve is given in parametric form by the function 
z(t) = (x{t),y{i)) as t varies, its digitization is essentially defined by the formula 

round z(t) = (round x(t), round y(t)) 

as t varies, where round(a) is the integer nearest a. 

We need to be careful, of course, when rounding values that are halfway between 
integers, because round (a) is undefined in such cases. Let us assume for convenience 
that the path z(t) does not go through any pixel centers; i.e., that z(t) is never equal to 
(m+ |, n+ 1) for integer m and n. (Exact hits on pixel centers occur with probability zero, 
so they can be ignored in the theorem we wish to prove. An infinitesimal shift of the path 
can be used to avoid pixel centers in general, therefore avoiding the ambiguities pointed out 
in Bresenham's interesting discussion [1]; but we need not deal with such complications.) 
Under this assumption, whenever we have x{t) = m + \ so that 'round is ambiguous, 
the value of round y(t) = n will be unambiguous, and we can include the entire line 
segment from (m, n) to (m + 1, n) in the digitized path. Similarly, when t reaches a value 
such that round x(t) = m but round y(t) = n or n + 1, we include the entire segment from 
(m,n) to (m,n+ 1). This convention defines the desired digitized path, round z(t). 

When the path z(t) returns to its starting point or begins and ends at infinity, with- 
out intersecting itself, it defines a region in the plane. The corresponding digitized path, 
round z(t), also defines a region; and this digitized region turns out to have a simple charac- 
terization, when we apply standard mathematical conventions about "winding numbers" : 
The pixel with corners at (m, n), (m + l,n), (m, n + 1), (m + 1, n + 1) belongs to the 
digitized region defined by round z(t) if and only if its center point (m + |, n+ |) belongs 
to the undigitized region defined by z(t). (This beautiful property of digital curves is fairly 
easy to verify in simple cases, but a rigorous proof is difficult because it relies ultimately 
on things like the Jordan Curve Theorem. The necessary details appear in an appendix to 
John Hobby's thesis [2], Theorem A.4.1.) 

Now we are ready to begin proving the desired result. The region defined by an angle 
at (xo, yo) with lines of slopes a/b and c/d can be characterized by the inequalities 

a(x - x ) - b(y - y ) > ; c(x - x ) - d(y - y ) > . 

(We may need to reverse the signs, depending on which of the four regions defined by two 
lines through (xo,yo) are assumed to be defined by the given angular path; this can be 
done by changing (a, b) to (—a, —6) and/or (c, d) to (— c, — d).) This region contains the 
pixel with lower left corner (m, n) if and only if 

a{m + \ - x ) - b{n + \ - y ) > ; c(m + \ - x ) - d(n + \ - y ) > . 

We can simplify the notation by combining several constants, letting a = a(xo — |) — 
Kvo ~ \) and 13 = c(x - \) - d(y - \): 

am — bn> a; cm — dn > (3 . 
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Let R(a, (3) be the digitized region consisting of all integer pairs (m, n) satisfying this 
condition; these are the pixels in the digitized angle corresponding to (xQ,yo). 

As noted above, it is safe to assume that the pixel centers do not exactly touch the 
lines forming the angle; thus we are free to stipulate that am — bn^a and cm — dn 7^ (3 
for all pairs of integers (m, n). However, if equality does occur, we might as well define 
the digitized region R(a, (3) by the general inequalities am — bn > a and cm — dn > (3, as 
stated, instead of treating this circumstance as a special case. Notice that R(a, (3) is equal 
to R{ \(3~\)', therefore we can assume that a and (3 are integers in the following discussion. 

Another corner point (x' ,y' Q ) will lead to parameters (ck',/3') defining another region 
R(a',(3') in the same way. The two regions R(a,(3) and R(a',(3') have the same shape if 
and only if one is a translation of the other; i.e., R(a, (3) = R(a',(3') if and only if there 
exist integers (k, I) such that 

(m, n) G _R(a, (3) (m — k,n — I) G R(a', (3') . 

Our main goal is to prove that the number of distinct region shapes, according to this 
notion of equivalence, is exactly \ad — bc\. 

Lemma. Let a, (3, a', f3' he integers. Then R(a, (3) = R(a',[3') with respect to slopes 
a/b and c/d if and only if a — a' = ka — lb and (3 — (3' = kc — Id for some integers (k, I). 

Proof. Assume that R(a,f3) = R(a',f3') with respect to a/b and c/d, and let (k,l) be the 
corresponding translation amounts. Thus we have 

J am — bn > a 1 J a(m — k) — b{n — I) > a' 1 

\cm- dn> 13) \ c(m - k) - d(n -l)>(3') 

for all integer pairs (m, n). Let a" = a' + ka — lb and (3" = (3' + kc — Id, so that 

J am — bn > a 1 J am — bn> a" 1 

\ c, -dn >(3 J \cm-dn> 13" J 

for all integers (m, n). This implies that a = a" and (3 = (3" . For if, say, we had a < a" ', 
we could find integers m and n such that am — bn = a and cm — dn > (3, because a and b are 
relatively prime; this would satisfy the inequalities on the left but not on the right. (More 
precisely, we could use Euclid's algorithm to find integers a' and V such that aa! — bb' = l. 
Then the values (m, n) = (aa' + bx, ab' + ax) would satisfy the left inequalities but not 
the right, for infinitely many integers x, because ad — be 7^ 0.) 

Thus R(ot, f3) = R(a' 7 (3') implies that a — a' = ka — lb and (3 — (3' = kc — Id. The 
converse is trivial. □ 

Let k and / be integers such that a = ka — lb. The lemma tells us that R(a, (3) = 
-R(0, (3 — kc + ld); hence every digitized region R(a, (3) has the same shape as some digitized 
region R(a', (3') in which a' = 0. 

It remains to count the inequivalent regions -R(0, (3) when (3 is an integer. According 
to the lemma we have i?(0, (3) = _R(0, (3') if and only if there exist integers (k, I) with 
= ka — lb and (3 — (3' = kc — Id. But ka = lb if and only if k = bx and / = ax for some 
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integer x; hence the condition reduces to (3 — (3' = bxc — axd = x(bc — ad). In other words, 
i?(0, (3) = i?(0, /?') if and only if (3 — (3' is a multiple of ad — be. The number of inequivalent 
regions is therefore |a<i — 6c|, as claimed. 

To complete the proof of the theorem, we must also verify that each of the equivalence 
classes is equally likely to be the class of the digitized angular region, when the intersection 
point (xo,yo) is chosen at random in the plane. The notational change from (xo,yo) to 
(a, (3) maps equal areas into equal areas; so we want to prove that the equivalence class of 
R(a, (3) is uniformly distributed among the \ad — bc\ possibilities, when the real numbers 
(a, (3) are chosen at random. Choosing real numbers (a, (3) at random leads to uniformly 
distributed pairs of integers ([a], \/3]). And if \a] has any fixed value and \(3~\ runs 
through all integers, the equivalence class of R(a, (3) runs cyclically through all \ad — bc\ 
possibilities. 

Q.E.D. □ 

A close inspection of this proof shows that we can give explicit formulas for the sets 
of intersection points (xo,yo) that produce equivalent shapes. Let D = \ad — bc\ and let 
Rj denote the shape corresponding to region i?(0, j) in the proof, where < j < D. Then 
the digitized angle will have shape Rj if and only if (xo, yo) lies in the parallelogram whose 
corners are (\ — bj/D,^ — aj /D) plus 

((b-d)/D,(a-c)/D), (-d/D,-c/D), (b/D,a/D), (0,0), 

or in a parallelogram obtained by shifting this one by an integer amount (m, n). 

In the special case a/b = 2/1 and c/d = —3/1, the shapes Rj are what we called Pj 
above; in the special case a/b = 3/(— 1) and c/d = —1/2, the Rj are what we called Qj. The 
shapes that appear in the digitized angles depend on the values of (xo mod 1, yo mod 1) in 
the unit square, according to the following diagrams: 




Notice that the parallelograms "wrap around" modulo 1, each taking up an area of 1/5. 

When the slope of either line forming an angle is irrational, the number of possible 
shapes is infinite (indeed, uncountable). But we can still study such digitizations by 
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investigating the shape only in the immediate neighborhood of the intersection point; 
after all, those pixels are the most critical for human perception. For example, exercises 
24.7-9 of The METRFONTbook [3] discuss the proper way to adjust the vertices of an 
equilateral triangle so that it will digitize well. 

The moral of this story, assuming that stories ought to have a moral, is probably 
this: If you want to bisect an angle in such a way that both halves of the bisected angle 
are visually equivalent, then the line of bisection should be such that reflections about 
this line always map pixels into pixels. Thus, the bisecting line should be horizontal or 
vertical or at a 45° diagonal, and it should pass through pixel corners and/or pixel centers. 
Furthermore, your line-rendering algorithm should produce symmetrical results about the 
line of reflection (see [1]). 

This subject is clearly ripe for a good deal of further investigation. 

Acknowledgments. I wish to thank the referees and the editor for their comments. In 
particular, one of the referees suggested the present proof of the theorem; my original 
version was much more complicated. 
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